addpath('\\128.112.86.75\lithium\Publications\Charge-Diffusion\DQMC');
addpath('\\128.112.86.75\lithium\Publications\Charge-Diffusion\high_temperature_expansion');

%%

% DateCell = {2017 08 21};
% Folders = [8,9,10,15,11,12,13,16,17,18,7];
% Times = [0,50,100,150,200,250,300,350,400,500,10e3];
% % Folders = [8,9,10,15,11,12,13,7];
% % Times = [0,50,100,150,200,250,300,10e3];
% ExcludePics = {[6,11],[11],[],[],[11,12],[],[37,39],[],[],[],[]};
% PeriodGuess = 5;

% DateCell = {2017 08 21};
% Folders = [19,21:27,20];
% Times = [0:50:300,350,10e3];
% % Folders = [8,9,10,15,11,12,13,7];
% % Times = [0,50,100,150,200,250,300,10e3];
% ExcludePics = {[],[],[],[],[],[],[],[],[]};
% PeriodGuess = 5;

%%%U/t = 8, HoldT = 2000ms
% DateCell = {2017 08 29};
% Folders = [28,29];
% Times = [0,10e3];
% ExcludePics = {[],[],[],[],[],[],[],[],[],[],[],[],[],[],[],[]};
% PeriodGuess = 5;

dset = DiffuseSet(DateCell,Folders,Times,PeriodGuess,ExcludePics,[],[],OnesFolder,ThreesFolder,ExcludeOnes,ExcludeThrees);
[FitParams,SErrs] = dset.fitSinglesVsDens(0,'DQMC',[1,1]);
dset.Temp = FitParams(1);
dset.TempUnc = SErrs(1);
% dset = DiffuseSet(DateCell,Folders,Times,PeriodGuess,ExcludePics,[],[],OnesFolder,ThreesFolder);
% dset.saveDat;